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Abstract 

The addition of an A4 family symmetry and extended Higgs sector to the stan- 
dard model can generate the tribimaximal mixing pattern for leptons, assuming 
the correct vacuum expectation value alignment of the Higgs scalars. Deviating 
this alignment affects the predictions for the neutrino oscillation and neutrino 
mass observables. An attempt is made to classify the plethora of models in the 
literature, with respect to the chosen A4 particle assignments. Of these models, 
two particularly popular examples have been analyzed for deviations from tribi- 
maximal mixing by perturbing the vacuum expectation value alignments. The 
effect of perturbations on the mixing angle observables is studied. However, it is 
only investigation of the mass-related observables (the effective mass for neutrino- 
less double beta decay and the sum of masses from cosmology) that can lead to the 
exclusion of particular models by constraints from future data, which indicates the 
importance of neutrino mass in disentangling models. The models have also been 
tested for fine-tuning of the parameters. Furthermore, a well-known seesaw model 
is generalized to include additional scalars, which transform as representations of 
A4 not included in the original model. 
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1 Introduction 



The experimental evidence of neutrino oscillations implies massive neutrinos, which 
contradicts the predictions of the standard model (SM). There are currently many ex- 
periments focused on precise measurements of the neutrino mass and mixing parameters: 
neutrino physics can be said to have entered the "precision era" . 

Global fits to the latest neutrino oscillation data (l]^ show that the leptonic mixing, 
or Pontecorvo-Maki-Nakagawa-Sakata, matrix f/pMNS is very close to the tribimaximal 
mixing (TBM) matrix 
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first proposed in Ref. [5]. Since the allowed deviations from TBM can only be small (not 
more than 10-15%), this mixing pattern represents at least a zeroth order approximation 
to lepton mixing [6] . It is completely different from the mixing in the quark sector, and 
has motivated extensive research into models of family symmetries [7|. 

Some of the discrete family symmetries used in the literature are^: ^4, S3, S4, T', 
A(27) and S(81); there are also models that employ continuous symmetries such as 
SU{3) or 5*0(3). The A4 models have a very economical structure in terms of group 
representations and field content. The most general mass matrix leading to TBM can 
be shown to be invariant under one of the group generators j7| (see the Appendix). 
Furthermore, the use of ^4 can be geometrically motivated: it is the symmetry group of 
the regular tetrahedron, and the angle between two faces is 26'tbm, where sin^ ^tbm = |- 
These characteristics have led many authors to construct and/or study models based on 
A4. Some models generate neutrino masses via effective dimension-5 operators, some 
apply the type I seesaw mechanism, whereas others use the type II, or the type I -|- 
II seesaw mechanisms. Table [l] is an attempt to classify the vast number of models,^ 
according to the chosen A4 assignment of the lepton doublets, lepton singlets and, if 
appropriate, the seesaw particles.^ The majority fall into the first four categories. 

Very often the TBM scheme is obtained only approximately, or with the cost of fine- 
tuning and/or various assumptions, such as vacuum expectation value (VEV) alignment. 
These alignments are chosen, or the models are explicitly constructed, in order to reach 
alignment, resulting in a certain mixing pattern (in this case TBM). However, cor- 
rections to the VEV alignment are expected, be it from renormalization, higher order 
operators, or the tree-level exchange of heavy fermions, for example. The aim of this pa- 
per is to study the effects of VEV-misalignment on the neutrino mass and lepton mixing 
observables. There already exist some numerical analyses (T2|[T3||3l]|54||59] focused on 
specific A4 models. In addition, the effects of higher order operators have been studied 



in A4 59 and S4 64 , 65 models, where the unperturbed VEV alignments predict exact 
TBM. This work emphasizes that observables related to neutrino mass (that is, the ef- 
fective mass for neutrino less double beta decay (Oz//3/3) and the sum of neutrino masses 
for cosmology) provide the best possibility to disentangle the models. Furthermore, and 
in contrast to previous studies, a more general VEV-misalignment is allowed for. The 



^ See the review in Ref. |7l for a list of references 
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^An earlier, much less complete classification can be found in 
^ There are also models that use the inverse and linear seesaw mec 



nanisms 
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type III seesaw mechanism |63|, with the same particle assignments as type D models 



, as well as the inverse 
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Table 1: Particle assignments of A4 models in the literature. Lepton doublets, charged 
lepton singlets and right-handed neutrinos are denoted by Lj, £^ and u^, respectively. A 
denotes the Higgs triplets in the type II seesaw mechanism. Models that also study the 
quark sector have the superscript *, those that embed ^4 into a grand unified theory 
group have the superscript *. 



Type Li A References 
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analysis in the present paper is focused on models of types A and B that predict TBM, 
as well as generalizations of these models to include more Higgs singlets. 

In this analysis,^ the chosen VEV alignment is modified by random complex devi- 
ations, perturbing the neutrino and charged lepton mass matrices from their original 
structure {M^, and Mi) to the perturbed ones, and M'^. The resulting neutrino mix- 
ing angles and mass-squared differences can be compared with current data (Table |2]). 
The well-known standard parameterization of the PMNS mixing matrix is 



PMNS 



C12C13 
— S12C23 ~ Ci2S23Sl3e 
'S12S23 ~ Ci2C23Si3e* 



S12C13 
C12C23 ~ ■5i2S23Si3e 
~Ci2S23 — Si2C23Si3e 




C23C13 






e*^3 







(2) 



''Other approaches to deviations from TBM can be found in Refs. 66 -73 
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Table 2: Best-fit values and allowed na ranges for the global three flavor neutrino 
oscillation parameters, from Ref. Isl. 



Parameter Ato^i (IQ-^ eV^) sin^ 6*12 sin^ 6*13 sin^ 023 | Ar7?4 1 (IQ-^ eV^) 

Best fit 7.67 0.312 0.016 0.466 2.39 

Icr range 7.48-7.83 0.294-0.331 0.006-0.026 0.408-0.539 2.31-2.50 

2cr range 7.31-8.01 0.278-0.352 <0.036 0.366-0.602 2.19-2.66 

3cr range 7.14-8.19 0.263-0.375 <0.046 0.331-0.644 2.06-2.81 



where Cij = cos^, Sij = sm9ij and ^12,6*13,^23 (0 < % < vr/2) are the three mixing 
angles. There are three phases in Eq. ([2]): S is the CP violating Dirac phase, and A2 
and A3 are Majorana phases, with < 6, A2, A3 < 27r. The two Majorana phases, A2 
and A3, do not affect the neutrino oscillation probability, but have an influence on the 
amplitude for Oz//3/3. 

One can also perform a "fine-tuning test" for each model, by examining the values 
that the mass matrix parameters must take in order to give the correct mass-squared 
differences, before perturbations are applied. Since these parameters generally originate 
from the product of some coupling constant with the VEV of a Higgs scalar, any close 
relationship between the parameters is highly unlikely, and could be evidence of fine- 



tuning in a particular model 13 



The paper is built up as follows: in Section [2] a type A model is introduced, it is 
examined for fine-tuning, the addition of Higgs singlets is discussed, and the model is 
analyzed for deviations from TBM; in Section [3] the same procedure is followed for a 
type B seesaw model. Section |4] presents the summary and conclusions, and for the sake 
of completeness there is a discussion of the group in the Appendix. 

2 The original Ma/Altarelli-Feruglio type A model 

In type A models, lepton doublets transform as 3, charged lepton singlets as 1,1', 1", 
and right-handed neutrinos are absent. In this case the neutrino mass usually comes 
from dimension-5 operators. Although Table [T] contains a long list of references for 
type A models, many of these works are phenomenological analyses of the same few 
models. The original model by Ma 19 is further developed in Ref. [s], where also an 



extra-dimensional solution to the vacuum alignment problem is provided. 
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The models in Refs. 19 and [8] employ the so-called Ma-Rajasekaran (M-R) basis 
for A4, in which neither M^, nor Mi is diagonal, but the product of the mixing matrices 
in each sector leads to TBM. In order to connect A4 models with the modular symmetry 



and thus the larger framework of string theory, the same model can be formulated 11 
in a different basis for A^ (the Altarelli-Feruglio (A-F) basis). In this basis the charged 
leptons immediately come out as diagonal, which means that the neutrino mass matrix 
is in the flavor basis, and is diagonalized by the TBM matrix. The two bases are simply 
related by a unitary transformation, and the multiplication rules differ (see the Appendix 
for details). 



^Note that the model in Ref. 19 contains 6 Higgs triplets, whereas the model in Ref. [s] uses 
dimension-5 operators. 
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Table 3: Particle assignments of the A-F A4 model. There is also an additional sym- 
metry, which decouples the charged lepton and neutrino sectors, and a f/(l) symmetry 
to generate the hierarchy of charged lepton masses. 



Lepton 




A, 
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2 


3 


e" 


1 


1 
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I" 


t" 
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Scalar 
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2 


1 




2 


1 




1 


3 




1 


3 




1 


1 



2.1 The original model in the A-F basis 

Along with the usual type A particle assignments for leptons (Table [T]), this model has 
two SM Higgs doublets, which are invariant under A4, as well as two triplets <y9 and 
Vj', and an singlet ^, all three of which are gauge singlets (Table |3]). These particle 
assignments, along with the A^ multiplication rules, lead to the Lagrangian 

[+ x,i\LL)" + x,i"{LL)'] + H.c. + . . . , (3) 



where (33) transforms as 1, (33)' transforms as 1', and (33)" transforms as 1", and y, 
Xa and Xd are dimensionless coupling constants. The notation in Eq. ^ follows the 



simplified description from Ref. 11 , where the Higgs doublet fields and h^, and the 
cut-off scale A are set to 1. Thus the term y^e^i^ipL) is in fact yee^{ipL)hd/ A, Xai{LL) is 
short for Xa^{LhuLhu) / A} and so on. The dots stand for higher dimensional operators 
- in this model these are suppressed by additional powers of the cut-off A, as long as 
the VEVs are sufficiently smaller than A. The two terms in parenthesis on the second 
line of Eq. (g come from additional Higgs singlets; these were not part of the original 



model, but one can show 13 that TBM can still be achieved with either two or three 



Higgs singlets in this model. This will be discussed in Sections |2.2| and |2.3 

Upon symmetry breaking, the VEVs of the Higgs singlet and triplets take the align- 
ments 

(0 = , {v) = {v. 0, 0) and {^') = {v\ v\ v') , (4) 
which lead to the charged lepton mass matrix 

M, = Vdj { 0' y, I , (5) 

where Vd is the VEV of the Higgs doublet hd- Thus the charged fermion masses are 

V V V 

rrie = yeVd-^ , nif, = y^Vdj^ , = yrVdj^ ■ (6) 
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Figure 1: Scatter plot showing allowed regions (for the 3a ranges of the oscillation 
parameters) in the a — d parameter space for the original A-F model, with one Higgs 
singlet, normal hierarchy and TBM. 



When only one Higgs singlet ~ 1) is present. 



mo 



/a+f 



V 



_d 

3 
2d 
3 



the neutrino mass matrix is 
-I \ 



a — 

2d 
3 



(7) 



with mo 



A ' 



2xa^ and d = 2xdj-, where Vu is the VEV of hu- The neutrino mass 



matrix is diagonalized by the transformation 



MM 



A 



^ diag(a + d, a, —a + d) 



(8) 



with U = f/xEM, as in Eq. Q. Thus TBM is achieved, and the neutrino masses are 
mi = mo(a + rf), m2 = moa and = mo(— a + d), which results in the sum-rule 
2m2 + ms = mi. Here the masses are understood to be complex, with the Majorana 
phases still attached. Note that with only one Higgs singlet it is impossible to get the 



inverted mass hierarchy in this model, as shown in Ref. 13 



It is interesting to note that in the case of one Higgs singlet, with the mass matrix 
in Eq. ([T]), some fine-tuning is required between the parameters a and d for the model 
to give the correct neutrino mass-squared differences |'13'|. This seems rather contrived, 
since a and d come from the products of different Yukawa couplings with the VEVs of 
the Higgs singlet ^ and triplet if', respectively. As can be seen in Fig. [T| if both a and 
d are real (as in Ref. [l3]), there is a linear relationship between the two parameters. 
If d is complex (as in this analysis), there is only a slightly greater allowed region in 
the a — d parameter space. Note that w.l.o.g., a can be chosen to be real. There are 
no perturbations applied in this case, and the parameter mo is set to 0.025 eV, the 
typical scale for the mass matrix of normally ordered neutrinos. In later cases, where 
the inverted mass ordering is studied, mo is fixed to 0.05 eV. The magnitudes of the 
parameters a and d (and later also c) are randomly varied in the range |a, c, d\ < 4, with 
their complex phases varying from zero to 2tt. 
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2.2 Two Higgs singlets 

Recall that only one Higgs singlet is introduced in the original model (Table |3|. However, 
in the framework of ^4 symmetry it is natural to take advantage of all representations 
of the group, and in this model it is also possible to achieve TBM with both two and/or 
three Higgs singlets 13 . In addition to the Higgs singlet ^, the singlets ^' and ^" can be 
introduced [Eq. (|3])], transforming as 1' and 1" under A4, respectively. The new singlets 
have the VEVs 

(0=Wc and (0=«6- (9) 

With only two Higgs singlets, there are three possible combinations (^, C,'] C,, C," and 
^"), but one can show 13 that only the singlets ^' and ^" can give rise to TBM. In 



this case, the resulting mass matrix is 



mo 



/2d 
' 3 



C + 



d 
3 

2d 
3 



3 » 



(10) 



where b = 2x 



and c 



2Xr 



Uc 

A • 



An additional condition for TBM is that b = c, which is 



a consequence of the necessary /i — r symmetry, and with this constraint the eigenvalues 
turn out to be mi = mo(— c + rf), 1712 = 2moC and m^, = mo(c+(i), with the new sum-rule 
m3 — nil = ^2- In this case, w.l.o.g., c can be chosen to be real. The scatter plots in 
Fig. [2] show that the c — d parameter space is quite tightly constrained (note that with 
additional Higgs singlets, the inverted mass hierarchy is now possible). 



2.4 



2.2 



1.6 



1.4 




-1.5 -1 -0.5 0.5 1 1.5 
c 




X d complex 
o d real 



0.5 



(a) Normal hierarchy (b) Inverted hierarchy 

Figure 2: Scatter plots of the c — d parameter space for the A-F model with two Higgs 



singlets [Eq. (10)], for normal and inverted hierarchy, with the condition b = c. In order 
to emphasize the difference between the complex and real case, the entire parameter 
space is not shown: in the complex case \d\ ranges up to 4 for the normal hierarchy and 
2.5 for the inverted hierarchy. 



^Although ji ~ T symmetry is required to get TBM, this forces one to impose the ad hoc relation 
b = c. 
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2.3 Three Higgs singlets 

If all three singlets C,' and C,") are present, the resulting mass matrix is 



mo 



/a+f 



v 



c + 



d 
3 

2d 
3 



3 > 

d 
3 

3 



(11) 



and the requirement for exact TBM is that a ^ b = c, which again reflects the necessary 
fi — T symmetry/ Here one can choose real a and complex c and d, w.l.o.g. This case 
is equivalent to the original Ma model in Ref. 



19 



and here there is more freedom 
in choosing parameters, as can be seen from the scatter plots of a — c — d parameter 
space in Fig. |3| There is basically no more tuning necessary in order to generate the 
correct mass-squared differences. The eigenvalues of the mass matrix in Eq. (fTTl), with 



a 7^ 6 = c, are mi = mo(ci — c + d), m2 = mo(a + 2c) and = mo(— a + c + d). 



2.4 Deviations from TBM in the A-F model 



The three mass matrices in Eqs. ([T]), ( [Io| ) and (11) are phenomenologically interesting, 
and will be numerically analyzed below. In order to study deviations from TBM, the 
VEV alignments of the Higgs triplets are perturbed, so that 



{^) = {v,vet,vet) and {^') = {v' ,v' {1 + e,),v' {1 + e^)) . 
Furthermore, in the cases of two and three Higgs singlets, 

& = c(l + e3) 



(12) 



(13) 



is defined in order to study the effect of changing the relative alignment of the Higgs 
singlets. Recall that the condition 6 = c is necessary for TBM in both the two and three 
singlet cases. 

With the above VEV-misalignment, the charged lepton mass matrix becomes 



Vd 



A 



(14) 



In the unperturbed case, the mass of each charged lepton la is ma = UaVdx [Eq- ([6])]. 
In this analysis, the mass scale Vdj^ is fixed to the tau mass, and each of the coefficients 
He, and Hr are varied randomly by 10% around their unperturbed values. Note that 
the charged lepton sector is unaffected by additional Higgs singlets, due to the presence 
of a Z3 symmetry (Table [s]) . 

The deviated neutrino mass matrix with one Higgs singlet is 



MP' 



/a+f 



mo 



-1(1 + 62) -f(l + ei)\ 



I -zy 3^ 

f(l + ei) 



(15) 



f(l + e2)/ 



^It can be shown 



13 that the conditions a 



c, a ^ b ^ c and a — c ^ b do not simultaneously 



give TBM and the correct mass spectrum. 
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(a) Normal hierarchy 




(b) Inverted hierarchy 

Figure 3: Scatter plots of the a — c — d parameter spaces for the A-F model with three 
Higgs singlets [Eq. (11)], for normal and inverted hierarchy, with the condition b = c. 



and with two Higgs singlets and is 



mo 



/f c(l + 63)-f(l + e2) 
c+f(l + ei) 



1(1 + \ 



c(l + e3) + f (1 + 62)7 



(16) 
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Figure 4: Scatter plots of (a) sin^ 6^23 against sin^ 6*13 and (b) Jqp against |f/e3| for tlie 
A-F model, normal liierarcliy, witli one (red circles), two (blue plus signs) and three 
(green crosses) Higgs singlets. The solid and dashed lines denote the la and 3a allowed 
regions, respectively. 



The most general case (three Higgs singlets) is 



mo 



/a+f 



v 



c(l + e3) 



c + 



2d 



-f(U 



e2j 



1 + ei) 



a — 



d 
3 

2d I 



(17) 



c(l + e3) + f (l + e2)y 



where the condition a ^ c still holds. One proceeds by diagonalizing the matrices in 



ch) 



Eqs. (15), (16) and (17). The perturbation parameters are in general complex, and the 



range [e^- '| < 0.3 is used throughout this work, with the phases varied freely. W.l.o.g., 
one can choose the parameters ei and e^*^ to be real. The other parameters a, h and d 
are varied as before, and mo is also fixed as described above. 

It is interesting to compare the deviations from TBM for different numbers of Higgs 



singlets, with the same perturbations applied to in each case [Eq. (14)]. Fig. |4] shows 
the results for the normal mass hierarchy (it is impossible to get the inverted hierarchy 
with one Higgs singlet). There are small differences, and in general one can conclude 
that with more singlets, greater deviation from TBM is possible. However, it is evident 
that if VEV alignment deviations are applied, the models deviate from TBM in a 
rather random fashion, and it is difficult to draw any firm conclusions from the plots of 
mixing angle observables. 

In contrast, the mass-dependent observables ^rrii (the sum of absolute neutrino 
masses) and (mee) (the effective mass for OuPP) allow for comparison between the three 
cases presented above (Fig. |5|, and can in principle be used to rule out some cases. 
These two observables are explicitly given as 



mj = mi + m2 + and {rriee) = \ U^i mi + U^^ "^2 + U^s "^3 1 



(18) 
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Figure 5: Scatter plots of {niee) against for the original A-F model, with one, two 

and three Higgs singlets. Black triangles (red squares) denote the normal hierarchy un- 
perturbed (perturbed) case; indigo circles (green crosses) denote the inverted hierarchy 
unperturbed (perturbed) case. 

It is useful to plot these two quantities against each other, in both the unperturbed 
and perturbed case. The solid black lines in Fig. [5] represent the allowed ranges for 
normal and inverted ordering, using the best-fit values of the oscillation parameters 
from Table [2| and varying the Majorana phases. The dotted and dashed lines include 
the 3a variation in the oscillation data, for normal and inverted ordering, respectively. 
The scatter plots display the results of the analysis discussed above. The deviations 
from TBM lead to more overlap between the normal and inverted hierarchies, with two 
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Higgs singlets. Increasing the number of Higgs singlets effectively increases the allowed 
range for both (mge) and ^ m,, and one can see that the three-singlet case corresponds 
to the most general TBM mass matrix [see Eq. (A19)]. To give one example of the 
consequences of Fig. |5| note from the middle left panel that if (mge) is experimentally 
determined to be less than about 10~^ eV, the case with two singlets and normal mass 
hierarchy can be ruled out. 

In general, i.e., without any model constraining the mass matrices, it is possible for 
(mee) to vanish for the normal mass hierarchy. In the case of one Higgs singlet, for 
example, vanishing (mee) means that the (1,1) entry of the mass matrix in Eq. ([T]) (or 
Eq. (15) in the perturbed case) is zero, i.e., a = 2d/3. Using the mass eigenvalues from 
Eq. (Sf^ it follows that the ratio of mass-squared differences is 



2a + d 
4a 



(19) 



which is inconsistent with the data (r should be close to 1/30). Perturbing the VEV 
alignment [Eq. (15)] and setting the (1,1) entry of Ml) to zero gives — g (l + e2 — 3ei), 
which can become sufficiently small. Indeed, in the plot (mee) can take values well below 
10~^ eV. Similar evaluations can be made for the other cases, in this and the next section. 



3 The Altarelli-Feruglio type B seesaw model 

According to the classification introduced in Table [T| type B models have lepton dou- 
blets transforming as 3, charged lepton singlets as 1, 1', 1", and right-handed neutrinos 
transforming as 3. Neutrino mass can be generated by the type I seesaw mechanism or, 
when weak scalar triplets are introduced, with the type I -|- II seesaw mechanism. 



3.1 The original A-F seesaw model 



The model in Section 2.1 can be extended by introducing right-handed neutrino fields 
u'^, transforming as 3 under 11 . The new Lagrangian contains all the terms in 
Eq. ([3]), along with the additional terms 

[+ xc7e'(^^V^)" + XBCil^'l^l'] + H.c. + . . . , (20) 



where y is a coupling constant.^ Most details of the model, including the VEV alignment 
in Eq. (|4]), remain the same, with the charged lepton mass matrix given by Eq. (|5|. The 
Dirac mass matrix is yVu times the identity matrix, and the Majorana mass matrix 
is 



R 



V 



_d 

3 
2d 
3 



3 » 



a — 

2d 
3 



A 



(21) 



^In Eq. (20 1 the compact notation of Eq. (|3| does not apply. 
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where a = 2xa^ and d = 2xD^- With the type I seesaw mechanism {M^ = {M^)'^Mj^^M^), 
the hght neutrino mass matrix is 



mo 



3a(a + d) 



f3a + d d 

2ad+(f 
d—a 

V ■ 



d \ 

d^-ad-So? 



d—a 
2ad+d^ 



(22) 



with mo = y As in Eq. ([8j), the matrix in Eq. (22) is diagonahzed by the TBM 
matrix, with eigenvalues mi = mo/(ci + d), 1112 = rriQ/a and = mQ/{—a + d), leading 
to the sum-rule 2/m2 + l/m^ = 1/mi. Note that the inverted hierarchy is possible in 
the seesaw version of this model, even with only one Higgs singlet. 

The fine-tuning test again shows that in order for the correct values of the mass- 
squared differences to be reproduced, the parameters a and d must take rather specific 
values, as shown in Fig. [6j There is a similar amount of tuning as in the A-F model 
without seesaw (see Fig. [T]). 



3.2 Two Higgs singlets in the seesaw model 



In the original A-F model, the addition of extra Higgs singlets still allows for TBM, with 



certain conditions (Sections 2.2 and 2.3). This idea can also be applied to the seesaw 



version of the model. Again, it is possible to introduce singlets ^' and transforming as 
1' and 1", respectively. However, just like the non-seesaw case, the singlet combinations 
C,, ^' and ^, ^" cannot give rise to TBM. That is only achieved with the two singlets ^' 
and resulting in the light neutrino mass matrix 

/-d^ -2{b + c)d-3bc ^b"^ + db + (c - d)d 3c^ + dc + {b - d)d\ 

3c^ - d^ - 2{b + c)d {c - d)d + b{d - 3c) , 
\ • • 3b'^ - d'^ - 2{b + c)d J 

(23) 



rriQ 
^(2) 




(a) Normal hierarchy (b) Inverted hierarchy 

Figure 6: Scatter plots of the a — d parameter space for the A-F seesaw model, with one 
Higgs singlet, for normal and inverted hierarchy. 
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d complex 
d real 




0.4 -0.2 
c 



0.2 0.4 



(a) Normal hierarchy 



(b) Inverted hierarchy 



Figure 7: Scatter plots of the c — d parameter space for the A-F seesaw model, with two 
Higgs singlet, for normal and inverted hierarchy. 



2xb^ and c 



2x 



The condition 6 = c is 



with m(2) = 3 (53 + c3 - (6 + c)d'^), b 
required for exact TBM, as before, and once again the Z3 symmetry means that the 
charged lepton sector is unaffected. The neutrino mass eigenvalues are mi = mo/(— c + 
d), m2 = mo /2c and = mo/(c + d), and the mass sum-rule l/ms — 1/mi = 1/1712 
applies. The scatter plots in Fig. [7| show the allowed regions in c — parameter space, 
and exhibit a similar level of tuning as the one singlet case (Fig. [6]). 

3.3 Three Higgs singlets in the seesaw model 

If there are three Higgs singlets present, the light neutrino mass matrix is given by 



(3) 



m 



(24) 



where the elements of the symmetric matrix M^^^ are 

M[f = 3a^ - _ sbc -2{a + b + c)d 



(3) 
12 

(3) 
13 

(3) 
22 



6 + c)rf + 3 [}? - ac) 



3c 



(3) 



M 



23 

(3) 
33 



{h + c- d)d + a{d - 3b) , 
3c^ -d^ - 3ab - 2{a + b + c)d , 
3a^ + da — 3bc + (6 + c — d)d , 
3b^ -d^ - 3ac - 2(a + 6 + c)d , 

(a + 6 + c)d'^). Once again, the condition a ^ b 



(25) 
(26) 
(27) 
(28) 
(29) 
(30) 



and m^'^^ = 3{a? — 3abc + b^ + c 
is required for exact TBM. In this case the neutrino mass eigenvalues become mi = 
^o/{cL — c + d), m2 = mo/ {a + 2c), = mQ/{—a + c + d), and there is more freedom 
in choosing parameters, as shown in the scatter plots of a — c — c? parameter space in 
Fig. [8} As in the non-seesaw model, for three singlets hardly any tuning is necessary. 
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(a) Normal hierarchy 




(b) Inverted hierarchy 

Figure 8: Scatter plots of the a — c — d parameter space for the A-F seesaw model, with 
three Higgs singlet, for normal and inverted hierarchy, with the condition b = c. 



3.4 Deviations from TBM in the A-F seesaw model 

The seesaw model can be analyzed for deviations from TBM due to VEV misalignment, 
following the procedure outlined in Section "lA above, with the same limits for the 
parameters. The VEV alignment is perturbed as in Eq. (12) and, for the cases of two 
or three singlets, as in Eq. (13). With the deviated Higgs triplet alignments of Eq. (12), 
the charged lepton mass matrix is again defined by Eq. (14), and the light neutrino mass 
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matrix is 

3mo 



m 



(1)' 



/ {{di - 3a)^ - 4(^24) -2dl + 3ad2 - did2 -2d| + Sada - d^d^, 

(cZ| - 6a(i3 - Ad-idz) ^a? + Sdio - 2d\ - ^2(^3 
y • • (o?| - 6ad2 - 4did2) 



m^^y = TJa^ - 9a {d\ + 2c?20?3) + 2 (d? - 3(^ic?2(^3 + 0^2 + c^s) > (31) 

with dx = d, d2 = d{\ + Ci) and d^ = d{l + £2). The deviated neutrino mass matrix for 
two Higgs singlets is 

Ml^y = ^ M(^)' , (32) 

where the elements of the symmetric matrix M^'^^ are 

M^i^' = d\- (3c + 2^2) (36 + 24) , (33) 

= 9b^ + Sdsb - 2dl + 3cdi - did2 , (34) 

Mif = 9c^ + 3d2C - 2dl + Sbdi - did^ , (35) 

Mg^' = (d2 - 3c)^ - 2di(36 + 2^3) , (36) 

M^f = -2dj - {d2 - 3c) (4 - 36) , (37) 

M^f = (4 - 36)2 - 2(^^ + 2(i2) , (38) 

and 

m^^y = -6(ii(d2(36 + ^3) + 3cd3) + (^3 - 36)^(36 + 2^3) 

+ {d2 ' 3c)2(3c + 2^2) + 2dl ; (39) 

for three Higgs singlets, the deviated mass matrix is 

M: = ^M(3)', (40) 

777, W 



where the elements of M^^^ are 



M^f = {di - 3aY - (3c + 2^2) (36 + 2^) , (41) 

M[f' = 96^ + Sd^b - 2dl + {di - 3a) (3c - ^2) , (42) 

M^f = 9c^ + 3d2C - 2dl + 3bdi - d^d^ + 3a(d3 - 36) , (43) 

M^f = -(3a + 2di) (36 + 2^3) + (^2 - 3c)' , (44) 

M2^J' = 9a' + 3dia - 2dt{ + 36(^2 - 3c) + 3cd3 - ^2(^3 , (45) 

= -(3a + 2rfi)(3c + 2^2) + (4 - 36)' , (46) 



and 



m^^y = 27a^ - 9a (96c + dl + 2^2(^3) + 276^ - 96 {2did2 + dl) 

- 6dids{3c + d2) + {d2 - 3c)'(3c + 2^2) + 2dl + 2d^ , (47) 

with (ii, (^2 and ^3 as defined above. The effect of changing the relative singlet alignment 
is studied by setting 6 = c(l + €3) in both the two and three singlet cases. 
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10"" 
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sin 6, 



0.4- 



0.3 




10 



One singlet 
Two singlets 
One singlet 



10" 



10" 



10" 



10" 



sin 6, 



(a) Normal hierarchy 



(b) Inverted hierarchy 



Figure 9: Scatter plots of mixing angle observables for the A-F seesaw model 11 , with 
one, two and three Higgs singlets, for both the normal and inverted hierarchy. 





(a) Normal hierarchy 



(b) Inverted hierarchy 



Figure 10: Scatter plot of Jcp against \Ue3\ for the A-F seesaw model, with one, two 
and three Higgs singlets, for both the normal and inverted hierarchy. 



Figuresj9] and [TO] sh ow scatter plots for the mixing angles, from diagonalization of 
Eqs. (31 ), p2]) and (HOj). Again, there are small differences, and in general the deviations 
from TBM can become larger with increasing number of singlets. However, there is little 
discriminative power with regards to the number of singlets, and also with respect to 
the model treated in Section [U 

In spite of this, the mass dependent observables, plotted in Fig. 11, allow some con- 
clusions to be drawn. For instance, in the one singlet case there is a distinct separation 
of normal and inverted hierarchy, and the normal hierarchy case is very different to the 
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Normal Inverted 

1/^0 
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E m. (eV) 



Figure 11: Scatter plot of {rriee) against for the A-F seesaw model, with one, two 

and three Higgs singlets, for normal and inverted hierarchy. See the caption of Fig. [5] 
for an explanation of the symbols used. 



non-seesaw model (Fig. [s]). As another example, if the normal mass hierarchy is favored 
by experiment and {rriee) is measured to be 0.05 eV, the upper left panel of Fig. 11 
shows that the seesaw model with one singlet can be ruled out. 



4 Conclusion 

The present paper is a study of deviations from TBM due to VEV misalignment in 
models. After an attempt to classify the vast amount of literature according to the 
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representations under which the lepton doublets, lepton singlets and seesaw particles 
transform under A^, two particularly popular examples from classes A and B have been 
focused on. The models have been checked for tuning and then generalized, in the sense 
that extra singlets, transforming under representations of A4 that are not used in the 
original models, are added. In general, the more singlets that are introduced, the less 
tuning there is. The most general VEV misalignment is allowed for, and the conse- 
quences for the lepton mixing observables are studied. Since these quantities have little 
discriminative power, the focus is shifted to the observables related to neutrino mass. 
The scatter plots of (mee) — ^ parameter space are different in each model, and allow 
one to distinguish different models, even after deviation of the VEV alignment. This 
is an indication of the importance of neutrinoless double beta decay and cosmological 
mass determination in disentangling neutrino mass models. 
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Appendix 
A A4 tetrahedral symmetry 



The following is an outline of the A4 symmetry group |9(|27)|74| , upon which the models 
in this analysis are based. 



A.l Introduction to A4 

Ai^ is the alternating group of order 4, and is also the group of all even permutations of 
four objects, isomorphic to the group of rotational symmetries of the regular tetrahe- 
dron. It is a finite, non-Abelian subgroup of SO (3) [35] and SU(3). A^ has 12 elements, 
which can be divided into 4 conjugacy classes with membership 1, 3, 4 and 4. The 
dimensionality theorem implies that there are 4 irreducible representations with dimen- 
sion dj such that ^ - (i| = 12. The only solution is c^i = ^2 = = 1 and 1^4 = 3, and 
the representations are labeled as 1, 1', 1" and 3, which means that there are three one- 
dimensional representations and one three-dimensional representation. The character 
table of ^4 is shown in Table H with u = e*^'^/^ the cube root of unity.^ 



A. 2 Different bases for A4 

There are two bases for A4 commonly used in lepton family symmetry models: the 
Ma-Rajasekaran (M-R) basis and the Altarelli-Feruglio (A-F) basis. 

^Note that uj = e^^^/^ = -1/2 + %/3/2 satisfies = and 1 + w + = q. 
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Table 4: Character table of A4, where n represents the number of elements in each 
conjugacy class. 



Class 


n 






1" 




Ci 


1 


1 


1 


1 


3 


C2 


4 


1 









C3 


4 


1 




u 





C4 


3 


1 


1 


1 


-1 



A. 2.1 Ma-Rajasekaran basis 

yl4 can be generated by two basic permutations S and T, given by S" = (4321) and 
T = (2314), where the generic permutation (1,2,3,4) — )■ (ni, n2, ^3, 714) is denoted by 
(nin2n3n4). It follows that 

S'^ = T^ = {STf = 1 , (Al) 

which defines a "presentation" of the group. The one- dimensional unitary representa- 
tions are generated by 



1 
1' 
1" 



S=l T= 1 , 

S=l T = e*2-/3 = 00 , (A2) 
S = l T = e*^"/^ = , 



and the three-dimensional unitary representation (in this basis) is built up from the 
generators 














1 






-1 













■ 


Vo 





-i) 




U 





0/ 



^= 0-1 , r= 1. (A3) 
\0 -1/ \1 0/ 

The 3x3 matrices of the natural three-dimensional representation 3 are: 
C2 

C4 : I -1 I , I 1 I , I -1 I , (A4) 




where each matrix can be generated by S and T in Eq. (A3). It is evident that the 



characters of the 3 representation (the last column of Table |4]) are simply the traces of 
the matrices in each class. 
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The multiplication rules are given by 



1x1=1, 

1' X 1" = 1 , 
i" X 1' = 1 , 
1' X 1' = 1 " , 
1" X 1" = 1' , 

3x3 = 1 + 1' + !'' 



3(js + 3g 



(A5) 
(A6) 
(A7) 
(A8) 
(A9) 
(AlO) 



where 3^^ and 3^ are "asymmetric" and "symmetric" combinations respectively. If 3„ ^ 
(oi, a2, as) and 3^ ~ 62, ^s) are two triplets transforming by the matrices in Eq. (A4) 



then the three singlets and two triplets in the product in Eq. (AlO) are 



1 = aibi + 02^2 + «3&3 , 
1' = aibi + UJ^a2h2 + ^0363 
1" = aihi + UJa2h2 + oJ^a^ih^ 



3i 



A. 2. 2 Altarelli-Feruglio basis 



(a2&3,«3&l,«l&2) 

(a362,ai&3,a2&i) 



(All) 
(A12) 
(A13) 
(A14) 
(A15) 



In the M-R basis, the generator S in Eq. (A3) is diagonal. However, one can also 



represent ^4 in a basis where T is diagonal, obtained through the unitary transformation: 



(A16) 














r 


= v^TV = 


(i 


UJ 



















ii 


-1 


2 2 


S' 


= V^SV = 


2 


-1 2 








2 


2 -1 



(A17) 



where 



V 




1 


1 


1 


1 






1 







(A18) 



It is known that the most general mass matrix leading to TBM, 

'A B 

1{A + B + D) 



m 



,TBM 




(A19) 



is invariant with respect to 5": (5")"^m™^S" = m™^. Note that the matrix V is the 
so-called "magic matrix" , which appears in some A4 models as the unitary matrix that 
diagonalizes the charged lepton mass matrix. In the 5", T' basis, the multiplication rules 



are identical to those in Eqs. (A5) - (AlO), but the product of two triplets gives the 
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composition of the following irreducible representations: 



1 = aibi + 0263 + ash , 
1' = 0363 + 0162 + a2&i , 
1" = a2&2 + aib-i + 0361 , 

3^ ~ ^ (2ai6i - 02^3 - 03^2 , 20363 - 0162 - 0261 , 20262 - 0163 - 0361) , 

^as ~ ^ ('^263 - 0362, 0162 - 0261, 0163 - 0361) . 



(A20) 
(A21) 
(A22) 

(A23) 
(A24) 



A. 3 Equivalence of the two bases 



The model presented in Section |2.1| can be formulated in the M-R basis, using the same 
particle assignments and the Lagrangian in Eq. (|3|). With the product decomposition 



rules in Eqs. (All) - (A15), and the triplet VEV alignment 

{ip) = {v,v,v) and {^') = {v',0,0) 
the charged lepton and neutrino mass matrices are 



(A25) 



Mf, 



V 



Ve 



Ve 



Ve 



.yr Vr^ Z/rW^ 




(A26) 



In this case, Mc is diagonalized by the magic matrix [Eq. (A18)], and is diagonalized 
by 






1 





1 





1 




1^ 









V2 


V2 



(A27) 



which combines with V in Eq. ( A18[ ) to give U^bm- The neutrino mass matrix in 
Eq. (A26) is equivalent to that in Eq. ([t]), with the change of basis induced by V. Thus 
the two bases lead to equivalent models, with the triplet VEV alignments in the charged 
lepton and neutrino sectors effectively swapped [compare Eqs. Q and (A25)]. Note that 
the change of basis will change the relative phases of the eigenvalues of M^. 
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